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Abstract 

Networks or webs of domain walls are admitted in Abelian or non-Abelian gauge theory 
coupled to fundamental Higgs fields with complex masses. We examine the dynamics of the 
^ ' domain wall loops by using the moduli approximation and find a phase rotation induces a 

- repulsive force which can be understood as a Noether charge of Q-solitons. Non-Abelian 

gauge theory allows different types of loops which can be deformed to each other by changing 
a modulus. This admits the moduli geometry like a sandglass made by gluing the tips of 
the two cigar-(cone-)like metrics of a single triangle loop. We conclude that the sizes of all 
loops tend to grow for a late time in general models with complex Higgs masses, while the 
sizes are stabilized at some values once triplet masses are introduced for the Higgs fields. 
We also show that the stationary motion on the moduli space of the domain wall webs 
represents 1/4 BPS Q-webs of walls. 



e-mail addresses: miiioru@df.unipi.it; f uj imori ,nagashi(§tli.phys . titech. ac . jp; 
nittaOphys-h . keio . ac . jp ; K . OhashiOdamtp . cajn .ac.uk; nsEikaiQth . phys .titech .ac.jp 



1 Introduction 



In various area of physics, many kinds of topological defects are expected to be produced at 
a phase transition via the Kibble mechanism [1]. More than two extended objects like cosmic 
strings or domain walls intersect or meet with angles in general, and therefore such a production 
inevitably results in networks or webs of these objects |2]. In the condensed matter physics 
several examples have been observed while it is not the case of particle physics, astrophysics or 
cosmology. Previously cosmic string junctions were suggested to be a seed of galaxy formation. 
Although such a possibility has been denied by a recent cosmic microwave background data, 
it is argued that they may still play a certain role. A domain wall network was proposed to 
explain dark matter /energy [3]. Future observation of such defect networks in our Universe 
certainly deserves to be explored. Usually dynamics of these networks have been studied by 
computer simulation. On the other hand, in the case of particle-like solitons such as monopoles, 
the analytic method of the moduli space (geodesic) approximation has been developed [H [5]. 
By this their low energy dynamics can be described as geodesies of the moduli space of these 
solitons. Therefore the determination of the moduli space is crucial for this task. In a previous 
paper [6] we have successfully constructed the moduli space of domain wall networks in a certain 
model which allows a supersymmetric generalization. Supersymmetry is expected to exist in the 
early Universe so this situation is realistic. 

In this paper we will work out the dynamics of domain wall networks using the moduli 
space approximation. We find that the sizes of all loops tend to grow for a late time in general 
models with complex Higgs masses, while the sizes are stabilized at some values once triplet 
masses are introduced for the Higgs fields. To the best of our knowledge this is the first example 
to discuss the dynamics of a composite system of solitons analyticallyll Our model here is a 
U{Nc) gauge theory coupled to Np Higgs fields in the fundamental representation, which can be 
extended to possess M = 2 supersymmetry. This model has been recently studied extensively 
because it allows many kinds of Bogomol'nyi-Prasad-Sommerfield (BPS) solitons, see [3 [HI E] 
for a review. Vacua are isolated and disconnected in theories with the number of flavors more 



^ We have analyzed analytically compressed walls which appear as limiting configurations of multiple parallel 
walls compressed to each otherlT^. Although these configurations may be regarded as composite solitons, they 
can be obtained as a smooth limiting point within a moduli space of multiple parallel walls, and do not present 
qualitatively new features unlike our present case of the 1/4 BPS webs of walls. 
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than the number of color, Np > Nc, and with non-degenerate masses for the Higgs fields [TO] . 
Parallel multiple domain wall solutions exist as 1/2 BPS states when the Higgs masses are real 
and non-degenerate. By introducing the method of the moduli matrix [HI [TT], analytic solutions 
of these domain walls were constructed in strong gauge coupling limit [12] (see [I3] for domain 
walls in U{1) gauge theory). This method was then applied to construct vortex solutions [H], 
vortex-strings stretched between parallel domain walls [15] , instantons inside a vortex-sheet [16] . 
Finally the most general analytic solutions of 1/4 BPS networks (webs) of domain walls have 
been constructed in models with complex non-degenerate Higgs masses [m m [i9]@ These 
solutions contain full moduli of a network with arbitrary numbers of loops and external legs of 
walls. Effective Kahler potential of 1/2 BPS solitons was constructed in the superfield formalism 
[22] and then it has been generalized to the case of domain wall networks [6] . The zero modes of 
external legs are non-normalizable and have to be fixed to discuss dynamics, while zero modes 
corresponding to loop size and associated internal phase are normalizable and appear as massless 
fields in the effective theory. We have constructed the effective Kahler potential and the metric 
of the simplest triangle loop in f/(l) gauge theory coupled with Np = 4 Higgs scalars and have 
found that the metric has a geometry between a cone and a cigar [6]. This metric is rather 
non-trivial since it is regular on the tip although it corresponds to shrinking loop. Therefore it 
is expected to describe a smooth bounce of the loop. 

In this paper we discuss the dynamics of loops of domain walls for 1) a triangle single loop 

in the simplest model of A'^c = 1 and A'^p = 4, 2) a double loop in the model with Nq = 1 and 

A^F = 6 and 3) a non-Abelian loop in the model with A'^c = 2 and A^f = 4. This paper is organized 

as follows. In section 2 we summarize the previous results on the construction of domain wall 

networks and the effective action on them. In section 3 we first investigate the dynamics of the 

single triangle loop. The moduli metric allows the f/(l) isometry whose orbit is parametrized 

by a Nambu-Goldstone mode of the flavor symmetry spontaneously broken by the configuration. 

Associated with this isometry, a conserved charge Q exists in the general motion of the moduli 

space. When Q = 0, a. motion of a shrinking loop is bounced and the phase is rotated with 

the angle vr after the loop completely shrinks. When Q 7^ a shrinking loop is bounced at the 

^ In A/" = 1 supcrsymmetric field theories, junctions of domain walls were previously found to be 1/4 BPS 
states preserving only a quarter of supcrsymmetry |20] . Exact solutions of wall junctions were constructed in |21j . 
See [6] for more complete references of domain wall junctions in A/" = 1 supersymmetric field theories. 
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minimum size of the loop determined by Q. In section 4 we investigate the dynamics of double 
loop. In this case there exist two conserved charges Q\ and Q2 correspond to the phases of the 
two loops. Both loops will grow after their sizes bounce at the minimum irrespective of their 
Q-charges. In section 5 we work out the dynamics of non-Abelian loop. In this case there exist 
two different configurations of non-planar webs with a non-Abelian loop, which can be deformed 
to each other by changing a modulus. The moduli space geometry looks like a sandglass which is 
made by gluing the tips of the two metrics of a single triangle loop. Each region of the sandglass 
metric corresponds to the configuration of each non-Abelian loop. Depending on the value of 
the conserved charge Q one configuration can or cannot change to the other configuration. In 
section 6 we turn on the triplet masses for the Higgs fields. In the context of field theory with 
eight supersymmetry charges this is possible in three space-time dimensions. We find the third 
masses induce the attractive force for the loop size whereas the Q-charge induces the repulsive 
force. Then the size of the loop is stabilized at some value where two kinds of forces are balanced. 
This mechanism to stabilize the size is the same as the one of the Q-lumps in nonlinear sigma 
models with a potential term [23]-[25]; the size of Q-lumps are stabilized by the Q-charge and 
the masses. Also, it was shown in [26] that 1/4 BPS dyon can be understood as stationary 
motion in the moduli space of BPS monopoles with a potential term induced by the masses. 
Dyonic instanton is also understood as stationary motion in the moduli space of instantons [27] . 
In the same way, our motion in the moduli space of the domain wall webs suggests BPS dyonic 
extension of domain wall webs. In fact it has been previously shown in [25] that the configuration 
of Q-domain wall webs is again 1/4 BPS (but not 1/8 BPS) and is stable. We reexamine this 
interpretation in this section. Section 7 is devoted to Conclusion and Discussion. Implication of 
our work to cosmology is briefly discussed. 

2 Effective Action of Domain Wall Networks 
2.1 BPS Equations for Domain Wall Networks 

Let us here briefly present our model (see [8] for a review), which admits 1/4 BPS webs of domain 
walls. We consider 3+1 dimensional M = 2 supersymmetric U (Nq) gauge theory with A'^p (> Nq) 
massive hypermultiplets in the fundamental representation. Here the bosonic components in the 
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vector multiplet are gauge fields Wm {M = 0, 1,2,3), the real scalar fields (a = 1,2) in the 
adjoint representation, and those in the hypermultiplet are the SU (2) doublets of the complex 
scalar fields {i = 1, 2), which we express as A^c x matrices. After eliminating the auxiliary 
fields, we obtain the bosonic part of the Lagrangian as 

2 



C = Tt 



- V, (2.1) 



V = 



■32 ■ 

a=l a=l ^ 

r,2 



(2.2) 



where we have defined = y (c^Iatq — {a'^Y ,-H^{H^y) with g the gauge coupling for U{Nc) 
gauge theory, and an SU{2)^ triplet of the Fayet-Iliopoulos (FI) parameters. In the following, 
we choose the FI parameters as c" = (0, 0, c > 0) by using SU {2)r rotation without loss of gener- 
ality. Here we use the space-time metric t/mn = diag(+l, —1, —1, —1) and Ma are real diagonal 
mass matrices. Mi = diag(mi, m2, ■ ■ ■ , 171^^), M2 = diag(ni, n2, ■ ■ ■ , n^p)- The covariant deriva- 
tives are defined as = c^a/S + i[WM, S], VmH"^ = {du + iWM)H\ and the field strength is 
defined as Fmn = -ipM, 'Dn] = duWN - OnVVm + i[WM, Wn]- 

If we turn off all the mass parameters, the vacuum manifold is the cotangent bundle over the 
complex Grassmannian T*GrNp,Nc [SO]- Once the mass parameters tua + iuA, (^ = 1, ■ ■ ■ Np) 
are turned on and chosen to be fully non-degenerate (m^ + in a 7^ ms + iriB for A ^ B), the 
almost all points of the vacuum manifold are lifted and only NpCnc = A^f!/[^cK-^f — ^c)'] 
discrete points on the base manifold Gri^fp^N^ are left to be the supersymmetric vacua [10]. This 
choice of the mass parameters breaks the SU{Np) flavor symmetry to U{1)'^^~^. Each vacuum is 
characterized by a set of Nq different indices {Ai, ■ ■ ■ , Aj^^) , 1 < < • ■ ■ < A]\r^ < Np, which 
we will often abbreviate as {Aj.) in the following. In these vacua, the vacuum expectation values 
are determined as 

(H'rA^ = ^cb^^A. {H^^'^) = 0, (S) = diag (mA, + zn^,, ■ ■ ■ ,mA^^ + m^^J , (2.3) 

where r is color index running from 1 to Nq, the flavor index A runs from 1 to N-p and E is the 
complex adjoint scalar deflned by S = Si + iTj2- 

The 1/4 BPS equations for webs of walls interpolating the discrete vacua (12. 3p can be obtained 
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by usual Bogomol'nyi completion of the energy density as [171 [18] 

Fi2 = 2[Si,S2], DiS2 = P2Sl, V^Ei+V2^2 = Y\ (2.4) 

ViH^ = H^Mi - V2H^ = H^M2 - T.2H^. (2.5) 

Here we consider static configurations which are independent of X3, so we set 9o = ^3 = and 
Wq = Ws = 0. Furthermore, we take = because it always vanishes for the 1/4 BPS 
solutions. The Bogomol'nyi energy bound is given by 

£>y + Zi + Z2 + daJa, (2.6) 

where the central (topological) charge densities which characterize the solutions are of the form 

3; = ^9„Tr(e°^S2P^Si), Z, = cdiTiEu Z2 = cd2TiT.2. (2.7) 

The topological charges are defined by 

(fx (Zi + Z2) , y= I d^x y. (2.8) 



Here corresponds to the energy of domain walls and Y corresponds to the energy of domain 
wall junctions. Since energy of domain walls means tension times the length of the walls, this 
quantity is divergent. On the other hand Y has a finite value, and we call this charge as the 
junction charge or the Hitchin charge. Note that the integration of the fourth term daJa = 
daTrlH^lMaH"^ — iJ^^S^)] in Eq. (12. 6 p does not contribute to the topological charges. 

The 1/4 BPS equations Eq. ([23]) and Eq. ([23]) [HI [H] can be solved as follows. Firstly, 
since the first two equations in Eq. (12.41) give an integrability condition for the two operators 
Va + So (a = 1, 2), Wa and can be written as 

S„ + iW^ = S-^d^S. (2.9) 

Here, S{x^,x^) G GL{Nq, C) is a matrix valued function. Secondly, Eq. (12. 5p can be solved as 

= S-'Hoe^''''+^'''\ (2.10) 

Here Hq, which we call "moduli matrix" , is an A^c x -^f constant complex matrix of rank Nq, and 
contains all the moduli parameters of solutions. Any sets of S and moduli matrix Hq related by 
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the following ^^-transformation are physically equivalent since they do not change the physical 
configuration: 

Ho -> VHo, S{x\x'^) VS{x\x'^), V e GL{Nc, C). (2.11) 

Finally, the last equation in Eq. (12. 4p can be converted, by using an A^c x -^c matrix valued 
function defined by 

n{x\x^) = SS\ (2.12) 

to the following equation: 

^{do.{n~'dM))= Inc - ^''^0, (2.13) 

where Qq = ^Hq e^^^^^^^^^^'^^^^ H^. This equation is called the master equation for webs of walls. 
Since HW — cl^^ = in vacuum regions, the solution r2(x^,x^) of the master equation should 
approach VLq near the vacuum regions. It determines S for a given moduli matrix Hq up to gauge 
transformations and then the physical fields can be obtained through Eq. (12. 9 p and Eq. (12.101) . 

There is a useful diagram to understand the structure of webs of walls, which is called the 
grid diagram [T71 [TB] . The grid diagram is a convex polygon in the complex plane Tr (S) (S = 
El + iT,2)- A vacuum point labeled by {Ai - ■ ■ A^^) correspond to the vertex of the convex 
polygon plotted at Tr (S) = Xl^i {^Ar + inAr)- ^ach edge connecting two vertices, there is 
a domain wall interpolating the two vacua and each triangle corresponds to a 3-pronged domain 
wall junction. Some examples of grid diagrams are shown in Fig.[I] and Fig.[2J In non-Abelian 




(a) triangle loop in Nq = 1, A'^p = 4 model (b) double loop in Nq = 1, A'^p = 6 model 
Fig. 1: Grid diagram and web diagram in Abelian gauge theory. 

gauge theory, two vacua with only one different label such as ( . . . A) and ( . . . B) can be 
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connected while two with ( . . . AB) and ( . . . CD) are forbidden to be connected. If there are 
several ways to connect the vacuum points, we obtain different configurations as shown in Fig.[2l 
By varying the moduli parameters, we can move one configuration to another one. 




(a) mass arrangement and vacuum points (b) non-Abelian loop in Nq = 2, Np = 4 model 
Fig. 2: Grid diagram and web diagram in non-Abelian gauge theory. 



One can easily read physical informations about domain walls and junctions from the grid 
diagram. The tension of the domain wall is proportional to the length of the corresponding edge 
of the grid diagram. More precisely, for a domain wall interpolating between vacuum {...A) and 
vacuum {...B), the tension is given by 

j^{...A){...B) =c\mA-mBl (2.14) 

where m^, friB are two component vectors such that friA = (m^,nA), ttib = ijnB.nB)- Fur- 
thermore the magnitude of the junction charge is proportional to the area of the corresponding 
triangle and its sign can be read off from the vacuum labels. If the junction interpolates three 
different vacua with labels such as {...A){...B){...C), this junction is called "Abelian junction" 
and its junction charge is given by 

y(...AK...BK...C) ^ _\^IABC][^ ^2.15) 

9 

where we have defined A[abc] as 

A[ABC] = ruA y< TUB + rfiB X rhc + rhc X rfiA, (2-16) 

which is twice the area of the triangle in the grid diagram. The junction charge above is negative, 
and can be interpreted as the binding energy of domain walls at the junction point. On the other 
hand, if the junction interpolates three vacua with labels such as {...AB){...BC){...AC), this 
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junction is called "non-Abelian junction" and its junction charge is given by 



lA 



Y{-Ab){...bc){...ca) ^ \-^[ABc\\ ^ ,2.17) 

9^ 

This is positive, and can be interpreted as the Hitchin charge of the Hitchin system. The details 
can be seen in 1181. 



In order to extract concrete informations from the moduli matrix Hq, it is useful to denote 
det i/^^'-^ = exp (a^^--) + ife^^""^), where ifQ^*"^ is an A^c x minor matrix whose elements are 
given by (//^^'■^)^* = {HoY^'. Defining the weight W^^''^ of the vacuum (Ar) = {A1A2 ■ ■ ■ A^^) 
by 

W^^^\x\x') = ^ {ruAy + riA^x^) + a^^^\ (2.18) 

r=l 

we can write the determinant of VLq as 

detfio = det(-Hoe'^''^^'-^''^'''^Hl \ = e""^'"''. (2.19) 

If only one of the weight W^"^''^ is non-zero, we can show that the configuration is the vacuum 
labeled by {A^). Since the solution of the master equation Vt is well-approximated by VLq near 
vacuum regions, we can estimate the position of the domain wall interpolating between vacuum 
{Ar) and vacuum {By) as a line on which the weights W^"^''^ and VV^'^''^ are comparable: 

Nc Nc 

y^{A.) _ y^{Br) ^ (^^^ _ x^ + Y^ {uA^ - nej x" + a^"^^^ - a^""^^ - 0. (2.20) 

Here the other weights should be sufficiently smaller than VV^"^''^ and W^^''^ . Hence the parameter 
Q(-4r> _ a^Br) |-]^g moduli matrix determines the position of the domain wall. Furthermore, one 
can see the angle of the domain wall is determined by the mass difference between the two 
vacua. Notice that the domain wall line Eq. fl2.20p is perpendicular to the corresponding edge 
of the grid diagram, see Fig.[Tl So the grid diagram gives us informations of the shape of the 
domain wall web as a dual diagram. A junction point at which three of domain walls get 
together can also be estimated by the condition of equating the weights of three related vacua 
as W^^'-^ ~ W^'^'-^ ~ W^'^'-^ 
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2.2 Effective Action of Domain Wall Networks 



Once we obtain the solutions of the BPS equations Eq. (12. 4p and Eq. (12. 5p . we can construct a low- 
energy effective theory on the world-volume of the domain wall networks. While all the massive 
modes on the background BPS solutions can be ignored at low-energies, moduli parameters (zero 
modes) as elements of the moduli matrix Hq can provide massless modes which will play a main 
role in the effective theory. Among these zero modes, we should promote only normalizable zero 
modes (pi to fields on the world-volume of the domain wall network as 

Ho{(l)')^Ho{f{xn), (2.21) 

where {fi = 0, 3) denotes the world- volume coordinates of the domain wall network. On 
the other hand, non-normalizable zero modes which change the boundary conditions at spatial 
infinities cannot be promoted to fields on the world-volume. 

In general, the master equation Eq. (I2.13P is difficult to solve. However we can obtain a general 
form of the effective Lagrangian for the moduli fields without solving the master equation, which 
have been constructed in [6]. It was found that the metric on the moduli space is a Kahler metric 
whose Kahler potential is given by 

1 



K(d>, 6) = / d'x 



clog det aoi(</', <P) + :-^IV(O-}(0, 0)a,aoi(</', 0))' 



(2.22) 



where Qsoi{4>, <P) is a solution of the master equation (I2.13p . In order to get this Kahler potential, 
one needs to solve the Gauss's law constraint for the world-volume elements of the gauge field 

V^F^o - I I?oS«] - 2y {H'VoH'^ - VoH'H'^) = 0. (2.23) 

We found [6] a generic form of the solution for the Gauss's law can be expressed by derivatives 
with respective to the moduli fields as 

w.ixn = I U.slM, Moxi^. 4>) - s;M MSsoM, 0)) , (2.24) 



where S'soi(0, 0) is given by Qsoi{<P,(p) = Ssoi{(p, 4>)Ssoi{4>, (pY and the variations are defined by 
Sfi = ^M0*^ ^^'^ ^ji ~ "^Ai^*^- ^^om this Kahler potential (I2.22p . the effective Lagrangian can 
be obtained as 

C'ff = didjK{(l), 0) dy'd^^^ = Kfj{<P, 0) d^<p'd^^K (2.25) 



The domain wall network in Fig.[T]-(a) (A'^c = 1, -^f = 4) has a single normalizable complex 
zero mode. To describe the zero mode, we can take the following moduli matrix without loss of 
generalitjc 



Ho 



% 1, 1, 



with 



r+ie 



(2.26) 



The complex parameter (p is the normalizable modulus parameter. One can easily see by looking 
at the weight of this system that its real part r changes the configuration of the triangle loop as 
shown in Fig.[3l and the imaginary part 9 corresponds to the phase of the loop. For sufficiently 
large r, the size of the loop is proportional to r. 




Fig. 3: The triangle loop configuration has four vacuum regions {A) [A = I,-- - ,4). We fix three 
complex moduli parameters which are related to positions of external walls. A unique normalizable 
mode is a zero mode which is related to the area of the region (4), namely the size of the triangle loop. 



The other zero modes, first three elements in Eq. (12.261) . are non-normalizable, and have to be 
fixed by boundary conditions when we construct the effective theory of the domain wall network. 
The effective Lagrangian for the size moduli have been already constructed in [6]. The Kahler 
potential and the metric are smooth everywhere in term of the complex coordinate 0, and the 
moduli space has a geometry between a cone and a cigar with a tip at = (r = — oo). It was 
found that the Kahler potential in the strong gauge coupling limit ^ oo is given as a sum of 
hypergeometric functions, see Eq. (3.30) of [6]. 

In particular the Kahler potential and the metric of the single triangle loop are given asymp- 

^ If we choose the V-transformation, the central position of the loop, and two relative phases carried by external 
walls, wc can always reduce the moduli matrix for the single triangle loop to the form (j2.26p . 
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totically for large |0| = by 
K 



4A 



[123] 



(log 101 ) T — \ + + (log 101 ) 



Qaia2Ciz g c \ ai a2 



4123] 



' - + 

aia20!3 g c \ as 





2 h 
- + - 


C V «3 




1^^231^ 1 


I'^Sl 




«2 



,(2.27) 



r _ 1 f\mi2\^ , |m23p , Irn^il'^ \ , , 2 , ,.2 



(c/r' + rfr), (2.28) 



where = — and ratios = ta^^ — ^abc^b x satisfying ai + 02 + «3 = 1- 

^^[123] 

The minus sign in Eqs. (12.271) . ( 2.281) is for the triangle loop in the U{1) gauge theory and plus 
sign in the f/(3) gauge theorjQ. The corrections to the asymptotic metric have been found to 
be exponentially suppressed. An interesting feature is that the above asymptotic metric can 
be understood as the kinetic energies: the first terms in the parentheses in Eqs. (12.271) . (12.281) 
represent the kinetic energies of domain walls and the second (with the =F sign in front) that 
of junctions. Since the lengths of domain walls are proportional to r, their masses and kinetic 
energies have linear dependence on r. On the other hand, the junction charges are localized at 
the junction points, and so their kinetic energies do not depend on r. This interpretation nicely 
explains the sign of the second term. See Eq. (I2.15P and Eq. (I2.17p . This result implies that the 
asymptotic metric for more complicated configurations can be also obtained by computing the 
kinetic energies of domain walls and junctions. We will often use this result in investigating the 
dynamics of various domain wall networks below. 

Before closing this section, let us briefly discuss another configuration closely related to the 
above one. When the vacua inside a loop are degenerate vacua, the loop acquires some internal 
modulilfl The asymptotic metric for the additional moduli exhibits another characteristic feature. 
Since the triangle loop requires at least 4 flavors that are non-degenerate, we need to take more 
than 4 flavors to examine a degenerate vacuum in the loop. Let us assume that mass parameters 
for external vacua are all non-degenerate 7^ rhs for A ^ B {A,B = 1,2,3,4), and the 
other mass parameters corresponding to the vacuum in the loop are all degenerate = rfiB 
for (4 < A, i? < A'^f). Such a loop configuration with the degenerate vacuum is described by the 



^ Here we take gauge coupling g and mass parameters for Higgs scalars (hypermultiplets) to be the same for 
two distinct U{1) and C/(3) gauge tlieorics with N-p = 4. They are dual in the sense that the number of vacua is 
equal and also their grid diagrams are congruent to each other. The duality becomes exact in the strong gauge 
coupling limit 00. 

^ Domain walls with degenerate masses were studied in [28l [29]. It was found that some Nambu-Goldstone 
modes for broken non-Abelian flavor symmetry are localized around (between) the domain walls and appear in 
the effective theory on them. 
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moduli matrix 



Ho = {1, 1, 1, cf>) , = (01, 02, ■ ■ ■ , 0iVp-3) . (2.29) 

In this case, there exist A^f ~ 3 complex normalizable zero modes: one of them corresponds to 
the size and phase of the loop and the others are zero modes associated with the vacuum moduli 
inside the loop. We can obtain the Kahler potential K in this case, if we replace |0p in Eg. (12 .2711 
by |0p = |0iP + ■ ■ ■ + |0Arp-3p. Therefore the knowledge of the Kahler potential for the Np = 4 
case gives a Kahler metric for this degenerate cases as 

= d^^d^M = 6,,K'{\<t)\') + 0,0,A'"(|0n, (2.30) 

where prime on K denotes differentiation with respect to See Appendix [Al for a concrete 

example. By differentiating the leading contribution (log |</>P)^ of this Kahler potential at asymp- 
totic region r = log 10 p ^ 1, we find that the Kahler metric contains not only terms proportional 
to r as in Eg. (12.281) . but also terms proportional to as shown in Eg. (]A.5|) . This feature shows 
that among moduli fields, there are massless modes with a support extending two-dimensionally 
over the entire vacuum region inside the loop in the web of walls as illustrated in Fig.Hl 




Xi Xi 

(a) size and phase moduli (b) vacuum moduli inside the loop 

Fig. 4: The densities of the metric for A'^c = ^,^F = 5 in strong coupling limit g ^ oo. At each 
point on the moduU space, the tangent space of the moduli space is orthogonally decomposed into the 
directions of size, phase and two vacuum moduli inside the loop. 
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3 Dynamics of Triangle Loop 



Since we have obtained the metric on the moduh space of the triangle loop, its dynamics can be 
discussed as geodesic motions on the moduli space. In order to avoid infinite volume of domain 
walls we may compactify the world- volume direction or simply dimensionally reduce the model to 
1+2 dimensions. Such model is obtained merely restricting the indices M, in the Lagrangian 
(12. ip to 0,1,2. The effective Lagrangian takes the form 




with w = r + iO m Eg. (12. 26 p . It is worth emphasizing that the moduli space is regular with 
positive curvature even when the loop shrinks completely. Fig.[S] shows the embedding of the 
moduli space into 3-dimensional Euclidean space. The moduli space has a U{1) isometry which 




Fig. 5: The moduli space of single triangle loop embedded in R^: The moduli space has a U{\) isometry 
which corresponds to the direction of the phase modulus. The other direction can be regarded as the 
direction of size modulus of the loop. The tip of the moduli space corresponds to the point (/> = where 
the loop shrinks completely. 

originates from a linear combination of three U{1) flavor symmetries. Correspondingly, there 
exists a conserved charge such that 



In terms of this conserved Q-charge, the effective Lagrangian can be rewrittero as 



We have performed a Legendre transformation of L with respect to 9. 
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Here the second term can be interpreted as a potential associated with the conserved charge Q. 
Note that the smoothness of the Kahler metric in terms of = e"' means the metric K^w is 
exponentially suppressed as Kyj^ oc e^*" -^^ for r = Re ^ — oo. We also know the asymptotic 
behavior Kyj^ oc r — > oo for r — > cxd. The typical form of the potential is shown in Fig.[6l 
We can find that the phase rotation produces the repulsive potential among the triangle loop. 




Fig. 6: Plot of the potential V = Q^/AK^y^ in — > oo limit (solid line) and asymptotic (r ^ 1) 
potential V = Q"^ /4,ar, a = 2a[i23] 010203 (hashed line) given in Eq. (j2.28p . The phase rotation produces 
the repulsive force among the triangle loop. 



This repulsive potential makes the loop to expand forever, namely the trajectory of the loop 
exhibits a runaway behavior. Although we are now considering the effective theory of domain 
wall networks, it is possible to consider domain wall networks with Q-charges in the original 
theory. The above runaway potential tells us that such configuration is unstable and no longer 
BPS. However, we will see in section [6l that a stable stationary point appears if we introduce 
another type of mass term (triplet mass). The corresponding configuration will turn out to be 
BPS , conserving a quarter of supercharges. 

Now let us return to geodesic motions on the moduli space. We introduce an integral of 
motion E as an integration constant as 

corresponding to the energy associated with the motion of the zero-modes r(t) and 9{t). By 
exploiting this conservation law of energy, we can obtain the solution of the equation of motion. 
The orbit of the geodesic for a given energy E is given by 

e-eQ = ± I dr , ^ (3.5) 
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and the time dependence of the size modulus r is given by 

t-t, = ± [ dr^Jh^k^=. (3.6) 

If we consider the motion in the direction of smaller values of r with Q = 0, the geodesic is a 
straight line in the complex (/)-plane and goes through the tip of the manifold = (r = — oo). 
This motion corresponds to the bounce of the loop, that is, after the loop shrinks completely, 
it tends to be larger with 180° phase rotation. In the case of Q 7^ 0, the repulsive force among 
the loop become stronger as the size become smaller, so that it prevents the loop from shrinking 
completely. Hence, there exists a minimum value of the size modulus r determined by 

^-''^ww\' ram J 

This implies that if the initial velocity of the size modulus ^ is negative, the loop shrinks to its 
minimum size r = r^^^ and then the velocity ^ changes its sign. 

The large size behavior can be investigated using the asymptotic metric Eq. f l2.28p . First 
note that the second term in Eq. (12.281) can be absorbed by shifting the parameter as r — 
r =1= aia2tt3 (|w^i2p/tt3 + |^23p/'^i + \^3i\'^/(^2) /g'^c- After the shift, the equation of motion for 



r can be solved as 



9-9, 







(3.8) 



'2aE 



t-to = ±^ {2aEr + Q') ^AaEr - Q\ (3.9) 
where a = ttx^ — . In the case of Q = 0, the above equation says r ~ ti. This reflects 

2A[i23] 010203 ) J 

the fact that the mass of the triangle loop is proportional to r and its velocity becomes smaller 
as the size of the loop becomes large. For the loop with Q ^ 0, the minimum size is given by 
'"min = Q"^ l^aE. The typical time dependence of the size modulus is shown in Fig.[71 Since 
our argument above is based on the asymptotic metric Eq. (I2.28P which is valid for r ^ 1, the 
solution Eq. (13.81) and Eq. (13. 9 p can be well trusted only when rmin = Q"^ /4,aE ^ 1. Of course, 
the energy E should be small enough so as not to excite the massive modes. 

The asymptotic potential V = Q'^/Aar can be interpreted as the shift of the energies associated 
with the walls composing the loop. This expectation can be confirmed by the following argument. 
For a domain wall with tension T, the rotation of its phase induces flavor charge density pq on 
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Fig. 7: The solution of the equation of motion for size modulus with Q ^ 0. The phase rotation 
produces the repulsive potential and the loop bounces back at rmin = Q"^ /^aE. 



the world-volume given by 



'T' - (f ) 

In addition, the phase rotation cause a shift of the tension as 



AT = T^l + pj/c^ - T ~ ^pj. (3,11) 

For the loop, each wall composing the loop becomes a domain wall with flavor charges by rotating 
the phase modulus 6. Since the tensions and lengths of the walls composing the loop are given 
by 

T^^'^) = c\mAl /(A4) = JM^^^, A = 1,2,3, (3.12) 

A[i23] aia2as 

the total flavor charge for the wall interpolating between (y4)-th and (4)-th vacua is 

Qa = l'""'"' = -T^ (3.13) 

From this expression, we find that the total charge Q = ^ Qa agrees with Eq. ( 13. 2p in large r 
limit and Qa are given by Qa = olaQ- Therefore the total shift of energy is given by 

A=\ A=l ^ ^ 

From the argument above, we can intuitively understand the reason why the phase rotation give 
rise to a repulsive potential as follows: Since the shifts of tensions of the walls are proportional to 
Pq and the lengths of the walls are proportional to r, the shift of the total energy is proportional 
to Pgr. This fact implies that if the charge density is constant, the total energy increases as the 
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size of the loop becomes larger. However since the conserved quantity is not the density pq but 
the total flavor charge Q, the density decreases in proportion to 1/r and hence the total energy 
decreases in proportion to 1/r as the size of the loop become larger. This is the reason why the 
phase rotation produces the decreasing repulsive potential. 

Next, let us consider the dynamics of a triangle loop with degenerate masses. The Kahler 
potential can be obtained by replacing |0p with \4)\^ = + ■ ■ ■ |</)Arp_3p in i^d^p) and the 
Kahler metric of the moduli space is given in Eq. (12. 30 p . All but the size moduli r = log|(/)| 
can be eliminated from the expression of energy by using the conserved charges associated with 
U{Np — 3) flavor symmetry. For example, the energy in the case of Np = 5 takes the form (See 
Appendix |A]) 

where |gp = Yll.=i QaQa and the conserved charges Q, qa{a = 1, 2, 3) are defined by 

« - (f - ' - - '^^^ (f - ^^V"'*') ■ 

Note that these conserved charges are related as Q = {(j)aa4>)qa/\(t>\'^ and satisfy an inequality 
< Igp. The second and third terms in Eq. (I3.15P can be interpreted as the effective potential 
V{r) associated with the conserved charges. For large r, this potential takes the form 

V-(r) = ^+ (3.17) 

The first term of the potential takes the same form as in the case of non-degenerate masses. 
Conversely, the second term is induced by the Noether charges associated with the vacuum 
moduli inside the loop. To understand intuitively the origin of the second term, we can use the 
same argument for the first term. What we should notice is only the fact that a part of the 
Noether charge which have no contribution to Q is supported by the two-dimensional vacuum 
inside the loop, while the charge Q has one- dimensional support on the walls. Therefore we 
can easily re-derive the behavior of the second term repulsive potential proportional to 1/r^. 
This potential can be also understood by exchange of the massless particles propagating the 
degenerate vacuum. 



17 



4 Dynamics of Double Loop 



We will consider the dynamics of double loop shown in Fig. [l]-(b) in this section. Unlike the pre- 
vious example, this configuration has two normalizable zero modes which are related to the sizes 
of the loops and their phases. By varying the sizes of the loops, we obtain various configurations 
of the loops. We will first explain the configurations of the domain wall web and then discuss 
the dynamics of the double loop. 

The model is U{1) gauge theory with six hypermultiplets and we choose six complex masses 
as follows (assuming a real positive value for m > 0): 

/ 3m ^/3m 3m ^/3m m m \ , , 

M = diag, ^— , --, , -, --j . (4.1) 

The solution of the BPS equations Eq. fl2.4p and Eq. 02.51) is characterized by Hq, which is now 
six component row vector 

Ho = ^(e«l+^^l,e''2+^''^e''3+^^^e''*+^^^e"5+^''^e°«+^^«). (4.2) 

Since the parameters and fej (i = 1, ■ ■ ■ , 4) are related to positions and phases of external walls, 
these will be non-normalizable modes if these are promoted to fields, and normalizable zero modes 
correspond to the parameters in the fifth and sixth components of Hq. They are related to sizes 
of the double loop and their phases. When we consider the effective theory of the domain wall 
network, we have to fix four complex moduli parameters and promote two normalizable modes 
to fields: 

Ho = y^(l, e3'"'/^ 1, e=^"^'/^ 0^(x^), ^\x^)) (4.3) 

with 0* = e^' = e^'^^^' {i = 1,2). For simplicity, we have chosen a somewhat symmetric set 
of four complex parameters for external walls. By varying the sizes of loops, we obtain seven 
different patterns of web configurations as shown in Fig. [8] and Fig. [91 Let us recall the vacuum 
assignment depicted in FiglT](b) which gives the grid diagram together with the web diagram A in 
FigJHl The fixed parameter / corresponds to the length between (1)(2)(4) junction and (2) (3) (4) 
junction, in other words, the length of (2) (4) wall (internal line) of configuration G in Fig.O We 
call the left loop surrounding vacuum (5) loop-1 and the right surrounding (6) loop-2. In the 
region A(B), both the loop-1 and the loop-2 appear as quadrangle loops (triangle loops). When 
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Fig. 8: Seven regions of moduli space correspond- Fig. 9: Configurations of double loop, 

ing to different patterns of configurations. 



the loop-l(2) grows and covers the junction (1)(2)(4) ((2) (3) (4)), the other loop-2(l) is eaten 
by the loop-l(2) as C(D) in Fig. [91 In the region E(F) the loop-2(l) vanishes and the triangle 
loop-l(2) exists. In the region G both the loop-1 and loop-2 disappear. 

Since the Kahler potential K is independent of the Kahler metric Kfj can be written as 
Kfj = -^-^K = ^-£j-^K. Then the effective Lagrangian takes the form 

L = KfJr^,r^)— — = KfJr^,r^) + , (4.4) 

^ dt dt ' \ dt dt dt dt J ' ^ ' 

where we have used Kf^ = Kf^. In this case, there exist two conserved charges defined by 
Qi = 2Kfj^. By using these conserved charges, the Lagrangian can be rewritten as 

L = Kf,{r\r^)^^ - \Knr'y)Q.Q„ (4.5) 

where K^^ is the inverse of the metric Kj:;. 

In the previous section, we have found the characteristic property of loop, that is, the loop is 
apt to become larger irrespective of Q = or Q 7^ 0. Therefore, we expect that the loops would 
become larger and sit in region A in Fig. [8] after sufficiently long time interval. Some examples 
of numerical solutions without flavor charges {Qi = Q2 = ^) are shown in Fig.[TOl For the 
initial velocities such that ^ ^ 2^ and ^ ^ |^ (Fig. [10]- (a)), the orbits of the solutions are 
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( a\ '^^^ < 9 dr^ dr^ > 1 dr^ (y^\ dr^ > q dr^ dr^ > dr^ ( \ dr^ ^ dr^ dr'^ > 1 dr^ 

\'^) dt ^ dt ' dt ~ 2 di y*^) dt dt -I dt r-u dt y^l dt dt ■• dt 2 dt 

Fig. 10: Numerical solutions of the equation of motion for double loop without flavor charges. The 
initial state has taken to be configuration A with the same loop size and some orbits for various initial 
velocities are shown in these figures. 

almost straight lines in r^-r^ plane and sit entirely in region A. For the initial velocities such that 
^ > 2^ and ^ > — ^ (Fig.[TU]-(b)), the orbits of the solutions first enter region C, namely 
one of the loops shrinks. Then, they bounce back at = — oo and return to region A. For the 
initial velocities such that ^ < ^"^^ IF ~ Ili" (Fig-[ini-(c)), one of the loops shrinks and 
bounces back at = — oo. Then, they enter region D and bounce at = — oo, namely the other 
loop shrinks to zero size. Finally, the orbits return to region A and the sizes continue to become 
larger. 

Next, let us consider the case of the double loop with the flavor charges. Fig.fTT] shows an 
example of the potential V = jK^^'QiQj and numerical solutions for the double loop with flavor 
charges. The potential increases rapidly outside region A and produces the repulsive force among 
the walls, so that any orbits of the solutions enter region A after sufficiently long time interval 
and continue to become larger. Other numerical simulations also demonstrate that region A is 
preferable. 

It is possible to know the asymptotic metric in region A for ^ ^ ml by computing the 
kinetic energy of domain walls. The kinetic energy of domain walls implies that the asymptotic 
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(a) Contour plot of the potential, (b) numerical solutions with flavor charges. 

Fig. 11: (a) Contour plot of the potential logV = log (^iK^*QiQj). (b) Orbits of numerical solutions 
with flavor charges. There exists a forbidden region where the potential energy exceeds the given total 
energy V{r) > E. 



metric ^ r"^ ^ ml is given by 



5v/3c 
2m2 



(8r^ - r'^)\dw^\^ + (Sr^ - r^)\dw^\'^ - (r^ + r^){dw^dw^ + dw'^dw^) 
Re /i+ + Re /i_ \djj-\^ 



(4.6) 



where /i± is defined by /i± = (1 ± A)u;^/2 + (1 =i= \)w'^/2, A = 3/v^. This form of the metric 
implies that the geodesic equation decomposes into two independent equations which can be 
solved as in the case of the triangle loop. 

FigHU appears to illustrate that the trajectory of the double loop configurations can bounce 
back at most only twice . However, this behavior is due to the particular mass assignment of the 
model, namely the center wall in configuration A is rather heavy. Let us consider smaller mass 
difference between two flavors corresponding to the vacua inside the loop, such as 

3m 



1 . I 3m . \/3m 
M = diag. I — , , 



\/3m m' m' 

— , —2 , — , I , m,' <^ m. 

2 2 2 2 ' 



(4.7) 



Then the mass of the center wall is much smaller than those of the other walls. Such mass 
assignment makes it possible that the double loop configuration bounces a lot of times. Moreover, 
if we consider the case of degenerate masses m' = 0, we can have configurations exhibiting as 
many "bounces" as one wishes. In this degenerate mass limit, the center wall is no longer visible. 
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rather, it spreads over the entire middle vacuum region in the loop. Correspondingly, the mode 
also spreads over the vacuum region as in the case of a triangle loop with degenerate masses 
(See Fig.H]) and it describes the degrees of freedom of the degenerate vacua inside the quadrangle 
loop. Furthermore, it is interesting to note that the infinitely many "bounces" in the degenerate 
mass limit naturally reduces to the repulsive force, which is described by a potential similar to 
the second term in Eq. fl3.17p with a different coefficient a. To describe the dynamics in the 
degenerate mass limit, we can use the expression for the energy Eq. fl3.15p . if the Kahler potential 
K is replaced by that for a quadrangle loop. 



5 Dynamics of Non-Abelian Loop 

We will next consider the dynamics of non-Abelian loop shown in Fig. [2] in this section. This 
configuration has four external walls, and also four internal walls which divide six vacua and 
constitute a quadrangle loop. After fixing the positions of external walls, one complex moduli 
parameter is left. The difference from the Abelian loop is that the moduli parameter controls the 
areas of two vacuum regions. First we will explain the configuration and moduli parameters, and 
then we will discuss the dynamics of the non-Abelian loop. For the details of the non-Abelian 
webs of walls, see [T8l. 



The model is U{2) gauge theory with A'^p = 4 hypermultiplets, and we choose four complex 
mass parameters as follows: 

T fm , 3m m . 3m\ , ^, 

Ml + iM2 = diag. \—-im, — , — + im, — — j . (5.1) 

The complex masses and vacuum points in the Tr(S) plane are shown in Fig. [2]-(a). The solutions 
of the BPS equations are characterized by 2 x 4 moduli matrix Hq. It is convenient to extract 
2x2 matrix if^^^^'^ defined by (iJ^^^^^^)^* = (ifo)'^% {s, t = l,2). Let us denote det i/^^^"^'^ as 

^{MA,) ^ exp(a<^i^2> + z6<^i^2>) ^ H^'^^'^'K (5.2) 
These parameters are not independent but satisfy the so-called Pliicker relation given by 

^(12>^(34> _ ^(13>^(24> ^ ^(14>^(23> ^ q_ (g^g) 

Each parameter a^^'^^^^ corresponds to the area of the vacuum region {A1A2) and 6^^i^2> 
the associated phase as before. In order to fix four external walls, we set four complex moduli 
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parameters as 



mL 



5(12} ^ ^(34> ^ ^(14> ^ ^(23> ^ Q_ 



(5.4) 
(5.5) 



The parameter L controls the positions of the external walls and the shape of the quadrangle 
loop as shown in Fig. [121 The remaining moduli parameters are r^^^^ and r^^^\ which determine 





x'-- 


--L/6 






<34> 
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L/2 






r x^^L/A 




' r/ <24> / 




= -L/A 






<12> 


<U> 

x'' = -L/6 


'l/3 






-l/3<s<l/3 1/3 <s 



(a) A configuration (b) Deformation of the loop in the web 

Fig. 12: The web diagrams in the parallelogram-type mass arrangement. Positions of the Abelian 

jmiction and the non-Abelian junction A'^* are given by Ai = {s — 1,— 1), = (|, I + s),^3 = 



(1^, 1), ^4 = (-1, -I - s), Ai = (1 + s, 1), A2 = (- 



2 

Off 



3' ■ 



|),Ar3 = (-4^,-l),A4 = (|,|-s) in unit 



,2 4 
^3' 3 

2 2 
3' 3 



the size of the loop. We introduce two complex parameters n, G C as 

a<i3>+,5(i3>^(^ + ,,)!!^^ (5.6) 
a<24>+,5(24>^(^_^)!!^. (5.7) 

The parameter u is fixed by the Pliicker relation Eq. fl5.3p . In the following, we take L sufficiently 
large, L ^ 1/m, so that the equation Eq. (15. Sp determines the parameter m as m ^ 1. Then the 
only parameter left is u, which we denote a.^ v = s + iO. The moduli parameter s controls the 
areas of two vacua (13) and (24), and three patterns of webs with a quadrangle loop appear 
as s changes, as shown in Fig.[T2l The parameter 9 is related to the Nambu-Goldstone mode 
corresponding to one of the broken flavor symmetries. 

Now let us discuss the dynamics of the non-Abelian loop. If we calculate the kinetic energies 
of domain walls separately in three configurations in Fig.[T21 the asymptotic metric on the moduli 
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space of the quadrangle loop is obtained as 





1 






1 



(5.8) 




In the outer two regions of the parameter s, the metric has linear dependence on s since the 



Fig. 13: (a) The metric evaluated numerically (solid line) and the asymptotic metric (jS.Sp computed 
from kinetic energies (dotted lines), (b) Embedding of the moduli space into the 3-dimensional Euclidean 
space. Here the metric is numerically evaluated in the limit 5 — > oo. The moduli space is non-singular 
since the curvature is finite everywhere. 

lengths of the internal walls depend linearly on s. We have observed the same feature in the case 
of the triangle loop in Eq. (12.281) . In the middle region, the linear dependence on s cancels out 
and the metric does not depend on s. Fig.[T3l-(a) shows the the numerically evaluated metric and 
Fig. [13]- (b) shows the shape of the moduli space isometrically embedded into the 3-dimensional 
Euclidean space. 

If we consider the motion without flavor charge from the outer regions of the parameter s to 
the direction of middle region, it goes through the middle regions and continue to go to the same 
direction. This motion corresponds to the motion of the loop changing the vacuum region inside 
the loop. If the conflguration has non-zero flavor charge Q 7^ 0, the potential term V = ^r^— 
will be induced in the effective Lagrangian. The typical form of the potential is shown in Fig. [Til 




(a) The metric of the moduli space (b) The moduli space embedded into R'^. 
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Fig. 14: The potential V 



If the energy E is greater than V^ax = ^7^^; the change of the vacua inside the loop can occur 
as in the case of Q = 0. However, if < V^ax, the quadrangle loop bounces back to be larger 
without changing the vacua inside the loop. 

Let us next compute the kinetic energies of junctions. The magnitude of the junction charge 
is proportional to the area of the corresponding triangle in the grid diagram in the complex Tr(S) 
plane. See equation Eq. fl2.15p and Eq. (12.1 7p . We show the areas of four junctions in Fig.[T51 The 




s < -1/3 -1/3 < s < 1/3 1/3 < s 

Fig. 15: The magnitude of the junction charge is proportional to the area of triangle in the complex 
Tr(S) plane which is dual to the junction point in the actual configuration. Abelian junctions Ai and 
non-Abelian junctions Ni (i = 1, ■ ■ ■ ,4) are illustrated in Figll2[ Non- Abelian junctions are denoted by 
shaded regions. 

sign of the Abelian junction is minus while that of the non-Abelian junction is plus. Although the 
total junction charge is zero in all three regions of the parameter s, the positions of junctions have 
different dependence on s. This causes the different velocities of junctions for a given value of ^, 
and the total kinetic energies of junctions can be non-zero. Since the Abelian junction transforms 
into the non-Abelian junction and vice versa at s = ±1/3, the kinetic energies of junctions are 
different in these three regions. These kinetic energies imply the additional contributions to the 
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asymptotic metric as 



3 f 1 

^^j' = ( 4 ) ^im (ds^ + dO^)^ ^ > 3' 

ds^ = 0, s^O, (5.9) 

6 Size Modulus Stabilization and Q Webs of Walls 

So far, we have seen that the sizes of the loops tend to become larger after sufficiently long 
time interval. In this section, we show that the sizes of the loops stabilize if the third mass 
parameters M3 are turned on in the Lagrangian (12.11) . One way to introduce M3 consistently 
with supersymmetry is the Scherk-Schwarz dimensional reduction from 1+3 dimensions to 1+2 
dimensions for the Lagrangian (12.11) . Then the third (twisted) mass parameter M3 is naturally 
introduced together with the third adjoint scalar E3 whose origin is the gauge field of the reduced 
dimension. 



6.1 Effective theory analysis 



Let us consider the triangle loop discussed in section [3] for simplicity. If we turn on a small third 
mass parameter such that M3 = diag (0, 0, 0, 7713), a potential Kn(r) is induced in the effective 
theory. We can show that this potential can be obtained from the 1+1 dimensional effective 



theory C = K^^) {d^rd'^r + d^Od'^O) by requiring r(t. 



X 



r{t), 6{t,x^) = 9{t) + m^x'^ and then 



reducing to 1-dimensional theory as 

2 



dr 
It 



fde 

\di 



Kn(r), 



Vm{r) 



(6.1) 



The asymptotic form of this potential obtained from the asymptotic Kahler metric (12.281) takes 
the form 



Vmir) = (msf 



(6.2) 



2A[i23] aia2a3 

which is valid for r ^ 1. This is a confining potential so that the loop shrinks and eventually 
shrinks to a point, if we do not turn on the fiavor Q-charges coming from the motion of the 
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phase. This potential can be interpreted shi 
mass parameter m^, the tensions of walls shift ai 



t of energies of the walls. Because of the small 



Ar<-^'« = ^(T(^.-i))2 + (m3c)2 - rf-*-"' « (6^3) 
Therefore the total shift of energy can be evaluated by using Eq. (13.121) as 

y AT<^'^>/<^'^> = (msf-r^ (6.4) 

^ 2A[i23] aia2«3 

When we turn on the flavor charge Q around the domain walls composing the loop given in 
Eq. (13. 2p . the effective Lagrangian can be rewritten as 



L = K^^{r)\^-j -VQir)-VUr), (6.5) 
VQ{r) + Kn(r) = + (m^fK^^ > \msQ\. (6.6) 

Let us remember that K^j^, vanishes in the limit of r ^ — oo and diverges in the limit of r ^ oo. 
Therefore, Vqi^r) = ^ increases as r ^ — oo while VrrSj) = {'^zf'^ww increases as r ^ oo 
asymptotically, and there is the minimum of the potential saturating the inequality in the last 
equation, as shown in Fig. [161 At the minimum, the value of the Kahler potential is related to 
the given Q-charge 



Q 



(6.7) 

ms 

Comparing this with the Q-charge for the unstable configuration in Eq. (13. 2p . we observe that 
the stable configuration has dO/dt = ±777,3. The relation (16.71) implies the size of the loop is 
stabilized at a certain value r = tq. For sufficiently small value of \m^\ and large value of the 
stabilized size tq takes a large value, so that it can be evaluated from the asymptotic potential 
as 

ro = -p^A[i23]aia2a3- (6.8) 
\m2\c 

Thus the third mass parameter M3 stabilizes the size moduli of loops by preventing the loops 
to expand forever. When the flavor charge Q becomes non zero, the configuration is stabilized 



^ The Schcrk-Schwarz dimensional reduction just introduces one more component of the energy density for the 
tension of domain walls in Eq. (|2.8|) . The tension of domain walls is formally still given by the same formula, 
proportional to the length of the mass vector: T^^'^'^ = c|m^ — ms|, except that the mass vector mA now becomes 
a three-vector in the three-dimensional grid diagram, after the Scherk-Schwarz dimensional reduction. 
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Fig. 16: The potential V for non-zero tliird mass parameter 7713 and flavor charge Q. The potential V 
is the sum of = {m^)"^ Kyj^j and Vq = 4^—- The potential takes the minimum |m3Q| at r = rg. 



at a finite size of the loop, rather than at the boundary of the moduli space corresponding to 
the complete shrinkage (vanishing size) of the loop. This stabilization mechanism is the same 
as the one of the Q-lumps in nonlinear sigma models with a potential term [23]-[25]; the size of 
Q-lumps are stabilized by the Q-charge and the masses. 

In the subsequent section we will show the BPS nature of this stabilized loop configuration 
from the viewpoint of original theory. There we will see the origin of the minimum of the potential 
(16. 6p which is always positive except for the case Q = or 7713 = 0. The value of the potential 
iQmsl at the minimum is shown to be equal to an increase in the BPS mass from the BPS loop 
with Q = ?Ti3 = 0. 

6.2 Q-domain wall web as 1/4 BPS soliton 

So far we have seen dynamics of Q charged domain walls and their networks mainly from the 
viewpoint of the low energy effective theory. Let's go back to the original theory and reanalyse 
the stable network with non-zero Q charges in more detail. It turns out that the configuration is 
a solution of another 1/4 BPS equations which are deformed from Eqs. (12.41) and (12.51) . In order 
to see it, let us consider supersymmetric model in c? = 2 + 1 with 8 supercharges in which there 
is an additional adjoint scalar E3 and third mass parameter M3 as mentioned above. 

It is convenient to write the mass matrix as M3 = mgH^ = diag (m\^rn\, ■ ■ ■ ,m^^), where 
we set Tr M3 = without loss of generality and (a = 1, 2, ■ ■ ■ , Np — 1) are the generators of 
U{1)^^~^, that is, the elements of the Cartan subalgebra of SU{Nf). The densities of conserved 
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charges of the U{1)^^ ^ symmetries are defined by 

In addition, it is convenient to define an electric charge density as 

pe = daTr (FoaSa) . 
Then the energy density can be written as 
£ = Tr 



(6.9) 



(6.10) 



Tr 



9 



fl' 9 \ ^ J 



+^ (F,o + ^^aSs)' + ^ (PoS, + z[S3, S,])' + 1 (P0S3)' 
5- 5- 5 



1\ |2 



+ y + Zi + Z2 + Ti&,pa + ^Pe + 5a 



> 3^ + Zi+Z2+m>„ + — pe + 5a^a 



where a stands for indices 1, 2 while a for 1, 2, 3, and we have used the Gauss's law 



(6.11) 



(6.12) 



The BPS equations are obtained by requiring the BPS bound to be saturated. Apart from the 
equations to determine time-dependence, we find the same five 1st order equations as Eqs. (12.41) 
and ([23D for the fields {Wa,^a,H^} 

F12 - ^[Si, S2] = 0, I?iS2 - V,^, = 0, I^iSi + I)2S2 = 4 {cIno - H^H^^) , 

ViH^ - HHh + J:iH^ = 0, V2H^ - HHl2 + J:2H^ = 0. (6.13) 

Then the solution, except for time dependence, can be written by using the solution of the 
master equation Eq. f l2.13p Q = SS"^ as Eq. (12. 9 p and Eq. (12.101) . Time dependence of all the 
fields including new variables {Wq, S3} are determined by additional equations [25] 



F„o + I^aS3 = 0, PoSa+i[S3,E„] =0, PoS3 = 0, 
VoH^ - iH^Ms + iJlsH^ = 0. 



(6.14) 



29 



If we choose a gauge such that Wq = —S3, these can be solved by replacing Hq as Hq Hq e*^^* 
and requiring the other fields to be independent of time 

iWi + El = S-^diS, iW2 + S2 = 3-^828, = S'^Ho e^'^^'+^^^^+'^^K (6.15) 

The spatial profile of the gauge field Wq (and adjoint scalar S3) are finally determined from the 
Gauss's law constraint Eq. fl6.12p . The solution takes the form of 

Wo = -S3 = -m^ (^£,S^S^~' + S~'^^S^ . (6.16) 

Here, a" (a = 1, 2, ■ ■ ■ , Np — 1) are the complex moduli parameters whose imaginary parts cor- 
respond to the Nambu-Goldstone modes of U{1)^^~^ symmetries, that is, a" appear in the 
expression of Hq as iJo(c") = Hq^c"" = 0)e'^"^". It is instructive to note a similarity between 
Eqs. f l6.16p and fl2.24p . Recall that the latter is the solution of the Gauss's law constraint for the 
configurations when we promote the moduli parameters of the stationary 1/4 BPS background 
to fields (functions of the world- volume coordinates including time t = x'^^^). Let us suppose 
that we make an Ansatz for the field cr" to depend only linearly on time as cT"(t) = a" + im'^t, 
namely the moduli matrix changes as ifg e'^"^" — > }J^Q^o.ia■'^+^mp) ^ This leads to dQa°'{t) = im'^ 
and then the solution fl2.24p corresponds to the solution (16.160 . 

For the solution Eqs. (16.150 and (16.160 . the Q-charges are determined by integrating the 
densities Qa = J <Px pa- Interestingly, the Q-charge which is a variable defined in the original 
theory can be directly related to the Kahler metric of the low energy effective theory (12.250 

Note that the Kahler potential K is independent of the imaginary parts of a". Since the right- 
hand side depends on the parameters contained in Hq., some of these parameters are fixed for given 
values of Qa- Therefore, those parameters are no longer moduli parameters and the configuration 
is stabilized. Especially, in the case of the triangle loop discussed in the previous subsection, we 
can show that the size parameter is fixed at the same value obtained as the minimum of the 
potential Eq. (16. 7p in the effective theory by taking, for example, a"" = ^0, 0, \J\'u^^ and mg = 
^0, 0, \J\f^ij with Ha=3 = ^^diag(— 1, —1, —1, 3) and the notation Q = ^J^Qa=z- Furthermore, 
the minimum value 1 7713(5 1 of the effective potential in Eq. (16. 6p precisely corresponds to the 
increment of energy bound mgQa in the last line of Eq. (16.110 . Note that since ^ at spatial 
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infinity in this case, the electric charge does not contribute to the total energy n Eq. (EUD tells 
us that the 1/4 BPS configuration requires both non-zero Q-charges and the third masses trig (or 
both of them to vanish simultaneously). If one of them is absent, the balance between them is 
lost and the configuration no longer is BPS. This is also consistent with what we found from the 
effective theory viewpoint in the previous sections. 

7 Conclusion and Discussion 

In this paper, we investigated dynamics of 1/4 BPS domain wall networks (or webs) in Abelian 
or non-Abelian gauge theories coupled with complex masses for Higgs fields in the fundamental 
representation. In the previous paper [6] we have obtained the effective action on the world- 
volume of the domain wall networks. In this paper we applied it to study the dynamics of the 
networks. Namely, We described the dynamics of the slowly moving networks as geodesies on 
their moduli space, namely with the moduli approximation. Only moduli parameters related to 
internal loops composed of several domain walls in the networks can be treated as massless fields 
in the effective Lagrangian. Other moduli are associated with the shift of external domain walls, 
which requires an infinite amount of energy and results in the change of boundary conditions. 

As concrete examples, we dealt with three different types of loops in Abelian or non-Abelian 

gauge theories. The first example is in Sec. [3] where the simplest configuration of the single 

triangle loop appears in the Abelian gauge theory with 4 massive Higgs fields. The metric of 

the moduli space for the single loop has a geometry between a cone and a cigar [6]. We found 

geodesies corresponding to any motion of a shrinking loop pass the tip (zero size of the loop). 

This means that the loop bounces back with vr rotation of the internal phase and eventually 

expands to an infinite loop. The second example is the dynamics of two loops in the Abelian 

gauge theory with 6 massive Higgs fields in Sec. HI There exist seven types of configurations shown 

in Fig. [9l We numerically showed that after sufficiently long time both of the two loops expand 

forever irrespective of the initial condition. As the two loops get larger, the system approaches 

to a system of two independent single loops. Our last example is the network including both 

Abelian and non-Abelian junctions which appears in the U{2) gauge theory with 4 massive 

* The Q-wall can be viewed as a capacitor with electric charge distributions on the two sides of the wall [JT] . 
Since these charges have opposite signs, the total electric charge vanishes. 
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Higgs fields [H], see Sec. [5l After fixing all the non-normalizable moduli, there remains only 
one complex moduli parameter s as the normalizable modulus which controls the areas of two 
vacuum regions. We found the metric of the moduli space whose geometry in Fig. [13] looks like 
a sandglass made by gluing the tips of the two metrics of a single triangle loop in Fig. [5l The 
geodesic of s is a one-way traffic from any initial value to an expanding loop with either s = +00 
or s = —00 (one or the other branches of the sandglass), depending on its initial velocity. Namely 
only one of the loop out of two loops remains after sufficient time. 

We also considered the dynamics of the web loop accompanied by a phase rotation in the 
internal direction, a U{1) isometry which originates from a linear combination of broken U{1) 
flavor symmetries. Conserved charges associated with the rotation are Q-charges of the domain 
walls composing the loops. These Q-charges give a runaway potential in the effective theory and 
exert a repulsive force between walls in the loop. Then the loops with Q-charges are generally 
unstable (non-BPS) and tend to expand forever. Thanks to the repulsive force, the geodesies 
bounce back before reaching the completely shrunk loops. The minimum sizes of the loops are 
determined by the given Q-charges. For the loops including both the Abelian and the non- 
Abelian junctions with the sandglass geometry in Fig. [13], the corresponding geodesic motion 
becomes bounce-back type or one-way traffic type depending on the total energy and the given 
Q-charge. 

By introducing the third masses for the Higgs fields, the effective theory of the loops acquires 
an attractive potential in contrast to the Q-charges. In the presence of the third masses, the loops 
tend to shrink. Coexistence of the Q-charges and the third masses stabilize the size of the loops. 
Then the size of the loop is fixed at some value where the attractive and the repulsive forces 
are balanced like the known stabilization mechanism of size moduli for the lumps due to the 
suitable potential accompanied by the Q-charge [23]-[25]. We also studied such configurations in 
the original theory, rather than in the effective theory on the world volume of the web of loops. 
Then we derived new 1/4 BPS equations which includes time derivatives and found a new BPS 
bound which is the sum of the topological charges of domain walls and their junctions and the 
additional masses coming from the Q-charges. General solutions of the 1/4 BPS equations with 
the Gauss's law are found. All the results found in the original theory are compatible with those 
found in the effective theory. 

Here we make several comments on possible extensions of the present work. 
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Global cosmic strings appear when global U{1) symmetry is spontaneously broken. It is 
well-known that there exists a repulsive force between two strings. Strings interact with the 
Nambu-Goldstone boson associated with the spontaneously broken U{1) symmetry, and the 
repulsive force was explained in terms of the Nambu-Goldstone bosons propagating in the bulk 
[2]. In the same way, it will be possible to explain the repulsive force induced inside a domain 
wall loop in terms of the Nambu-Goldsonte bosons. A feature different from the case of cosmic 
strings is that the Nambu-Goldstone modes in this case of a domain wall loop are normalizable 
and therefore appear in the low energy effective action of the loop as shown in this paper. 

An extension to a supertube [32] is interesting and may have some impact on string theory. 
With a Noether charge density, domain walls with arbitrary shape [33] were constructed as a field 
theory realization of a supertube. Our work should be extendible to a BPS supertube junction. 
Such a solution may suggest a junction of membrane tube in M-theory. 

It is interesting to explore applications of our results to cosmology. Our theory is supersym- 
metric and all stationary configurations discussed in this paper are BPS and are stable. In the 
early Universe supersymmetry is expected to be unbroken. Therefore our results can be applied 
when gauge and global symmetry are broken above the supersymmetry breaking scale. Vacuum 
regions inside domain wall loops are considered to be bubbles. Our results imply all bubbles 
grow for a late time in theory with complex masses if the size of the Universe is infinite. In 
this respect, it is worth generalizing our work to domain wall webs in a finite size space. In this 
case, zero modes of external legs of walls become normalizable and are promoted to fields in the 
effective theory. Dynamics of webs is not restricted to loops and will become richer. If we do 
not restrict ourselves to supersymmetric Universe, we can allow triplet masses of Higgs fields 
even in four space-time dimensions. In this case, the bubble (loop) sizes are stabilized. Growing 
bubbles with (without) a Q-charge will be stabilized (shrink) after supersymmetry is broken and 
the triplet masses are induced. 
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A Triangle loop with degenerate masses 

Let us consider a triangle loop for A'^c = 1? = 5 case in which the masses for fourth and fifth 
flavor components are degenerate m4 = m^. There exist four Killing vector fields C,o, (c^ = 
1, 2, 3) on the moduli space, which are given by 

eo = t(t>'di + (c.c), = t {(To),, <j>^d, + (c.c). (A.l) 

These Killing vectors originate from the U (2) flavor symmetry which rotates the fourth and fifth 
flavor components of H. Not all of them are independent. Instead they are related as 

^0 = ^ {M) ia. (A.2) 

The tangent space of the moduli space can be orthogonally decomposed into the direction of size of 
the loop tr- = (j)^di+{c.c.), phase of the loop tg = and two directions of the vacuum moduli inside 
the loop tj = c'}^a {I = 1,2). Here the coefficients d} are defined by [(j)a"-(f)) = 0, c^c} = Sjj. 
The norms of these vector fields are given by 

WUf = \\ter = 2|0p(/r(|c/>p) + \cf>\'K"{\cf>\')) = l^K, (A.3) 
lltif =||t2f = 2|0pA-(|0p) = |-K, (A.4) 

Here r = log \ 4>\ can be interpreted as the size of the loop. For large r, K — > 3^^^,^ 0,^020,3 ^^ 
the norms Eq. flA.3l) . flA.4p become 

= Wtef - -r^ = - TT^ (A.5) 

A[i23] aia2a3 ^[123] aia2a3 

These asymptotic form of the norms and their dependence on the size of the loop r show the 
fact that the metric densities for size and phase moduli have a one-dimensional support on the 
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edges of the loop, while those for the vacuum moduli have a two-dimensional support extended 
fully inside the loop (See Fig.H]). Since the moduli space has isometrics generated by the Killing 
vectors .^o and C,a, there are conserved Noether charges defined by 

Q = (0, ^o) = (0V' - 0V') , (A.6) 
qa = {<P,Q = tK,-(^}^{aa)\<P'-<P'^H^a)j), (A.7) 

where (p = (p^di + (c.c.) and ( , ) denotes the inner product with respect to the metric of the 
moduli space Kij. These conserved charges are related as 

Since the tangent vectors t,., tg, ti, t2 are orthogonal, the time derivative of the moduli parameters 
<p can be wirtten as 

^ = %Mi, + %Mt„ + %Mt, + <;^i, (A.9) 

II ''II ll^^'ll ll^lll ll^2|| 

Then the energy in the effecitve theory can be written as 

E = -(6 6) = 1 ( (0^ tr)' ^ (0, te)' ^ (0, h)' ^ (0, hY' 



2^ ' 2 \^ ||t,||2 ||t,p ||ti||2 \\t,\\^ 

-'- q2 TV- • 2 I / 1 1 \ /n2 I 9'a?a 



Here we have used (0, U) = r Ht^lp and c?c5 + + (0cr"0)(0crV)/l</>l^ = 5"^- 
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